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Abstract 

The effects of the final state interaction phenomenon known as colour reconnection are investigated 
at centre-of-mass energies in the range y/s — 189-209 GeV using the OPAL detector at LEP. Colour 
reconnection is expected to affect observables based on charged particles in hadronic decays of W + W~. 
Measurements of inclusive charged particle multiplicities, and of their angular distribution with respect 
to the four jet axes of the events, are used to test models of colour reconnection. The data are found 
to exclude extreme scenarios of the Sjostrand-Khoze Type I (SK-I) model and are compatible with 
other models, both with and without colour reconnection effects. In the context of the SK-I model, 
the best agreement with data is obtained for a reconnection probability of 37%. Assuming no colour 
reconnection, the charged particle multiplicity in hadronically decaying W bosons is measured to be 
(nS?) = 19.38 ± 0.05(stat.) ± 0.08(syst.). 
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1 Introduction 



Hadronic data in e + e~ collisions can be characterised by event shape distributions and inclusive 
observables such as charged particle multiplicities and momentum spectra. Measurement of the de- 
tailed properties of the hadronic sector of W + W~ decays allows the question of "colour reconnection" 
(CR) [1] to be addressed experimentally, in addition to providing tests of Monte Carlo models. The 
decay products of the two W boson decays have a significant space-time overlap as the separation 
of their decay vertices at LEP2 energies is small compared to characteristic hadronisation distance 
scales. In the fully hadronic channel this may lead to new types of final state interactions. Colour 
reconnection is the general name applied to the case where such final state interactions lead to colour 
exchange between the decay products of the two W bosons. A modification of the colour flow in this 
way could have a significant influence on the measured mass of the W boson, as first noted in [1]. It 
is therefore essential to ascertain whether or not such effects are present in W decays. As described 
in [2], a precedent is set for such effects in colour suppressed B meson decays, e.g. B — ► J/ipK, where 
there is "cross-talk" between the two original colour singlets, c+s and c+spectator. 

There is general consensus that observable effects of such interactions during the perturbative 
phase are expected to be small [2]. In contrast, significant interference in the hadronisation process 
is considered to be a real possibility. With the current understanding of non-perturbative QCD, such 
interference can be estimated only in the context of specific models [1-10]. Other final state effects 
such as Bose-Einstein correlations (BEC) between identical bosons from different W decays may also 
influence the observed event properties. 

This paper presents two different measurements which are sensitive to colour reconnection effects. 
The inclusive properties of W + W~ decay products have been measured [11-13] and found to have 
limited sensitivity to colour reconnection using the data available at LEP2. Characteristic observables 
such as the inclusive charged particle multiplicity in W + W~ — ► qqqq events, and its centre- 

of-mass energy dependence have been widely used to quantify the effect of colour reconnection in 
W + W~ events [1-6,8-10], and are therefore studied in this paper. As in [12], the hadronic part of 
W + W~ — > qq£T>£ events is compared with W + W~ — > qqqq events, while the leptonically decaying W 
is excluded. 

More recently, all LEP collaborations have concentrated on studies of "particle flow" [14-16], a 
generalisation of the well-known "string effect" [17] to the four-jet case of W + W~ — > qqqq, as models 
predict [18] this has a larger sensitivity to colour reconnection. This analysis compares the density of 
charged particles in two regions: the first, between pairs of hadronic jets originating from the same W 
boson, and the second, between pairs of hadronic jets which originate from different W bosons. In the 
absence of colour exchange between the two W bosons, the particle density is expected to be larger in 
the first region. Colour reconnection would lead to a migration of particles into the second region, in 
addition to a change in the total multiplicity. All data in the range 189-209 GeV are studied using the 
particle flow method in this paper, which supersedes previous OPAL analyses on the subject [11,12]. 

This paper is organised as follows: Section 2 summarises data and Monte Carlo models used, 
Section 3 describes the inclusive charged particle and particle flow analyses and Section 4 the estimation 
of systematic effects. Sections 5 and 6 discuss the results and draw conclusions. 
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2 Data Selection and Monte Carlo Models 



This paper uses data corresponding to an integrated luminosity of approximately 625 pb recorded 
during 1998-2000 with the OPAL detector, which is described fully elsewhere [19]. The data are 
separated into samples at six centre-of-mass energies, varying between approximately 189 GeV and 
209 GeV. Those accumulated above 202.5 GeV are considered at a single luminosity weighted mean 
centre-of-mass energy of 206.0 GeV. The selection criteria and distribution of data by W + W~ final 
state, luminosity and centre-of-mass energy are given in [20], with a total of 5401 W + W~ — > qqqq and 
2757 W + W™ — > qqfi^ candidates selected. Only W + W~ — > qq£Pi events in which the charged lepton 
was identified as an electron or a muon are used, to ensure that the charged particle multiplicity of 
the hadronically decaying W is well understood. 

It is essential to have reliable selection of charged particles in the detector in this analysis. Charged 
particles may have up to 159 hits in the principal tracking chamber, the jet chamber. Tracks used in 
the analysis are required to have a minimum of 40 hits in the | cos 9\ region 1 in which at least 80 are 
possible. At larger | cos0|, the number of hits is required to be at least 50% of the expected number 
and not fewer than 20, corresponding to a fiducial acceptance of | cos#| < 0.96. Tracks must have a 
momentum component in the plane perpendicular to the beam axis of greater than 0.15 GeV/c, and 
a measured momentum p of less than 100 GeV/c. For each track, the point of closest approach to 
the collision axis is found, and the distance between this point and the average interaction point is 
required to be less than 2 cm in the r-cj) plane and less than 25 cm in z. Clusters of energy in the 
electromagnetic calorimeter are required to have a measured energy greater than 0.10 GeV if they 
occur in the barrel region of the detector (| cos#| < 0.82), and greater than 0.25 GeV if they occur in 
the endcaps (0.82 < | cos6>| < 0.98). 

Most samples of Monte Carlo events used in this paper include detailed simulation of the OPAL 
detector [21] and of initial state photon radiation and have been passed through the same selection 
and analysis procedures applied to the data ("detector level"). A second class of samples does not 
include initial state photon radiation or simulation of the detector and allows all particles with lifetimes 
shorter than 3 x 10~ 10 s to decay ("hadron level"). 

Detector level samples were generated for a default set of physics processes at all centre-of-mass 
energies considered. Additional samples were generated for systematic studies as described in Sec- 
tion 2.4. Hadron level samples were produced for all variants of W + W~ events 2 considered (different 
hadronisation and colour reconnection models) at all centre-of-mass energies. The Monte Carlo event 
generators used to simulate the physics processes are described in the remainder of this Section, with 
emphasis on the CR models themselves. 

The effects of colour reconnection are implemented in several W + W~ event generators, and three 
groups of such models are studied, namely those of Sjostrand and Khoze (SK) [2,22], and those 
implemented in ARIADNE 4.11 [5,6] and in HERWIG 6.2 [7,8]. Events from all CR models used in this 
paper have been generated in conjunction with the electroweak generator KORALW 1.42 (KW) [23]. 
For the SK models, samples of events were generated such that they are identical to the conventional 
KORALW W + W~ — > qqqq events up to the end of the parton shower. This allows the construction of 
samples with an arbitrary fraction of reconnected events (including detector simulation) for the SK-I 
model, and also improves the statistical precision of studies using the SK models, such as estimation 

1 The OPAL coordinate system is defined such that the origin is at the geometric centre of the jet chamber, z is 
parallel to, and has positive sense along, the e~ beam direction, r is the coordinate normal to z, 9 is the polar angle 
with respect to +z and 4> is the azimuthal angle around z. 

2 In this paper, "W + W~ events" implies doubly-resonant W pair production diagrams, i.e. f-channel v c exchange and 
s-channel Z°/7 exchange, referred to as "CC03" in [14]. 
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of CR bias in measurements of the W boson mass [20]. For other CR models, the electroweak process 
was generated using KORALW and then a single set of events was hadronised by each of HERWIG 
and the ARIADNE models. The same KORALW events were also hadronised using the conventional 
QCD models of Section 2.4. All models considered have been tuned to describe Z° data, as described 
in [20]. 

2.1 SK CR Models 

The SK models are based upon the Lund string picture of colour confinement, in which a string is 
created that spans the decay product partons associated with each W. These strings expand from 
the respective decay vertices and subsequently fragment to hadrons. Before this occurs, at most one 
reconnection is allowed between sections of the two strings. The main scenarios considered are called 
type I (SK-I) and type II (SK-II) in analogy to the two types of superconducting vortices which 
could correspond to colour strings. In the SK-I scenario, the two colour flux tubes have a lateral 
extent comparable to hadronic dimensions. The probability for reconnection to occur is given by 
preco = (i _ exp(— Vki)), where V is the space-time integrated product of the maximum colour field 
strengths of the two overlapped W strings and kj is a free (dimensionless) strength parameter. In the 
SK-II scenario, the strings have infinitesimally small radii and a unit reconnection probability upon 
their first crossing. A third scenario considered, SK-II', is similar to SK-II but reconnection is only 
allowed to occur at the first string crossing which would reduce the total string length of the system. 

As described in [22], the only tuning necessary for these models is to ensure that the JETSET 
hadronisation model gives a good description of Z° data: parameters governing the behaviour of the 
reconnection model are not adjusted to fit data. Therefore, the same parameters were used as for 
the corresponding sample of non-reconnected e + e~ — > W + W~ events. The parton shower cut-off 
parameter, Qq, to which the predictions of the SK-II and SK-II' models in particular are sensitive, 
is set to 1.9 GeV in the OPAL tune [24] of the JETSET hadronisation model. The fractions of 
W + W~ —> qqqq events in which reconnection occurs at y/s = 199.5 GeV, P TCCO , are predicted to be 
17.2% for SK-II and 16.1% for SK-II'. As the fraction of events reconnected varies with kj in the SK-I 
model, two illustrative values of kj are given for comparison in figures and tables: ki = 0.9, giving a 
fraction of reconnected events, P ICCO ~ 34.3%, comparable to that used in [2], and an extreme case of 
ki = 100 (P rcco ~ 98%). The latter will be referred to hereafter as SK-I with 100% CR. 

Samples of these three models including simulation of the detector were generated at y/s = 
188.6 GeV, 199.5 GeV and 206.0 GeV. 

2.2 ARIADNE CR Models 

The second set of CR models is contained in the ARIADNE Monte Carlo program. They may be 
considered as extensions of the earlier partonic dipole model 3 [4], as both models were implemented 
using the ARIADNE Monte Carlo program and the same criterion is employed in the reconnection 
ansatz to determine whether reconnection is allowed. Perturbative QCD favours configurations with 
minimal string length in hadronic Z° decays [14]. When the partons of two W bosons are separating 
and strings are being formed between them, it is plausible that configurations corresponding to a 
reduced total string length are favoured. In the reconnection model of ARIADNE, the string length 
is defined in terms of the A measure, which may be viewed as the rapidity range along the string: 

3 In [4], at most one reconnection was allowed per event and possible reconnections between the decay products of a 
single W were not implemented. 
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A = ln(mf/mV), where rrn is the invariant mass of string segment i and m p sets a typical hadronic 
mass scale. Reconnections are only permitted if they satisfy the constraints of colour algebra and 
also lead to a reduction in the total A of the system. The first model is a variant of ARIADNE in 
which rearrangement of the colour flow is allowed but is restricted to the decay products of each W 
separately. This is referred to as AR-1. The second ARIADNE model, referred to herein as AR-2, is 
the same as AR-1 but in addition allows reconnections between the two W bosons for gluons having 
energies, E g < r\y, while the third ARIADNE model, AR-3, does not impose such a restriction. 
As gluons emitted with E g > Ty? ~ 2 GeV are perturbative in nature and have been shown to be 
radiated incoherently by the two initial colour dipoles [2], the model AR-3 is disfavoured on theoretical 
grounds. In addition, the ARIADNE colour reconnection models have been shown to be disfavoured 
by Z°/7 — ► qq data [25]. The AR-3 model is therefore not considered further in this paper. 

The way in which CR is implemented in AR-2 leads to an artificial difference relative to the AR-1 
model which is not directly related to reconnections between the two W systems. In AR-2, the dipole 
cascade (ordering in transverse momentum) is run in two stages from the maximum allowed gluon 
energy down to the cutoff: once down to E g = T\y allowing only reconnections within a single W 
system, and then a second time allowing E g < T\y and cross-talk between the two W systems. In 
AR-1, the dipole cascade is carried out in a single stage without any interruption down to the cutoff. As 
noted in [5], the AR-2 scheme is not strictly consistent with the assumptions of ordering in transverse 
momentum in the dipole cascade model and this leads to the observable differences between AR-2 and 
AR-1 referred to above. To ensure that differences between these two models are only due to inter-W 
reconnections, the dipole cascade in AR-1 is modified to run in two stages with an interruption at 
E g = T W [26]. 

As the same tuning of model parameters is used for both AR-1 and AR-2, and no colour flow takes 
place between the two W bosons in events in the AR-1 model, the AR-1 model serves as the no-CR 
model when estimating the expected sensitivity of the measurements to colour reconnection. The AR-1 
model is also used as an alternative model when estimating systematic effects in the hadronisation of 
the W decay products. 

Samples of the AR-1 and AR-2 models including simulation of the detector were generated at 
\fs = 199.5 GeV, with the fraction of W + W~ — > qqqq events in which reconnection occurs being 
approximately 49.4% for AR-2. 

2.3 HERWIG CR Model 

The third model is contained in the HERWIG program and provides an alternative CR model based 
on cluster hadronisation. In the cluster model, quarks and gluons from the perturbative parton shower 
evolution are combined locally into colour singlet objects called clusters which have (relative to strings) 
low mass and small space-time extent, each cluster decaying directly into a small number of hadrons. 
In the CR version of this model, an alternative pattern of cluster formation is implemented after 
the parton shower and gluon splitting phase. In this, new associations of partons into clusters are 
considered where they would lead to a smaller space-time extent of the clusters. When such viable 
alternative parton-cluster associations exist, they occur with a probability equal to | (= 1 /-^colours)- 
A sample of events including simulation of the detector was generated at y/s = 199.5 GeV, with 
reconnection occurring in approximately 23% of W + W~ — > qqqq events. 
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2.4 No-CR Monte Carlo Samples 



The models used are the same as those in [20], where a more detailed description can be found, and 
all are generated at the detector level. Samples of W + W~ events without colour reconnection effects 
are simulated at all centre-of-mass energies using the KORALW Monte Carlo, with fragmentation 
carried out using the JETSET 7.408 model. At three centre-of-mass energies (188.6 GeV, 199.5 GeV, 
206.0 GeV), additional samples are used in which the underlying W + W~ production process is sim- 
ulated by KORALW, while the fragmentation of a given set of four fermions is performed by either 
HERWIG, ARIADNE or an older parameter set of JETSET derived from tuning the model to OPAL 
inclusive event shape data, as used in [27]. Hereafter, these samples are referred to as JETSET, HER- 
WIG, ARIADNE and "old JETSET", respectively. Samples of W+W" events including Bose-Einstein 
correlations are simulated using the LUBOEI model [28] in PYTHIA 6.125 [3]. 

The dominant backgrounds in the W + W~ — > qqqq channel are e + e~ — > Z°/7 — > qq with ra- 
diation of energetic gluons, and four-fermion final states, primarily e + e~ — > (Z° /^)(Z° /j) — > qqqq. 
Backgrounds in the W + W~ — ► qqfi^ channel are significantly lower from all sources and receive small 
additional contributions from e + e~ — > Z°/7 — > r + r~, simulated using the KK2f version 4 Monte Carlo 
program [29]. For completeness, the small background represented by hadronic two-photon events is 
simulated using the PHOJET [30] and HERWIG event generators. 

Four samples of two-fermion processes e + e~ — ► Z°/7 — ► qq are simulated at each centre-of-mass 
energy, to allow systematic uncertainties to be estimated: in three samples the hard process is gen- 
erated using KK2f, with fragmentation of the quarks performed by each of JETSET, HERWIG and 
ARIADNE, while in a fourth sample both the hard process and fragmentation are generated by 
PYTHIA 6.125. 

Four-fermion processes are modelled using the KORALW 1.42 Monte Carlo, which contains matrix 
elements calculated by grc4f 2.0 [31]. The complete four-fermion samples are divided into two cate- 
gories: "WW-like four-fermion events" , corresponding to final states which could have been produced 
by diagrams involving at least one W boson, and "ZZ events" , which are the complementary sample, 
not all of which involve two Z bosons. 

In this and similar analyses, (Z°/7)(Z°/7) — > qqqq events are considered as background. In 
general, the susceptibility to the effects of CR in such events is expected to be comparable to that in 
W + W~ — > qqqq. A more complete treatment of such events would require implementation of the CR 
models in a full four-fermion generator, thereby obviating the need explicitly to subtract ZZ four quark 
final states. Although this is not done, it is expected to have only a small effect as the background 
level from (Z°/j)(Z°/j) — > qqqq events is low (less than 5% in the particle flow analysis). 

This four-fermion "background" is constructed from the difference between the predictions of two 
classes of events generated using the KORALW model: one containing the full set of interfering four- 
fermion diagrams (WW-like four-fermion events and ZZ events), the other containing only the W pair 
production diagrams. 

Alternative modelling of the four-fermion process with a more complete treatment of so-called 
0(a) photon radiation has been studied using samples generated by the KandY generator [32] which 
incorporates KORALW 1.51 [32] and YFSWW3 [33]. 
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3 Data Analysis and Correction Procedure 



The measurements of the inclusive charged particle multiplicity and of the particle flow are described 
below. The former is a fully inclusive measurement and data are corrected for the effects of finite 
detector resolution and acceptance, whereas the latter compares the predictions of a variety of models 
with the data at detector level. 

3.1 Inclusive Charged Particle Multiplicity 

The analysis of charged particle multiplicity follows the unfolding procedure described in [12]. The 
distributions of particle multiplicity and of the scaled charged particle momentum, x p = p/Ebeam, 
where -Ebeam is the beam energy, are used to measure the mean charged particle multiplicities in 
W + W~ — > qqqq events {(n^)) even ts, and in W + W~ — ► qq£u£ events {(n^")), and their difference 
(A(n ch ) = (nl q )-2(n^)). 

Figures 1 and 2 show the uncorrected multiplicity and x p distributions for W + W~ candidate events 
before background subtraction. The background predictions are the sum of all other Standard Model 
processes, as described by the models outlined in Section 2.4. The data are described reasonably 
well by all W + W~ models including conventional QCD processes alone, and by those including CR. 
Integration of the x p distribution is used for the principal measurement of mean charged particle mul- 
tiplicity as it has slightly lower estimated systematic effects than the direct multiplicity measurement, 
which is therefore used as a cross-check. 

The x p distribution is corrected for contamination using a bin-by-bin subtraction of the expected 
background, based on Monte Carlo estimates. Corrections are then applied for finite acceptance and 
the effects of detector resolution, using two samples of e + e~ — > W + W~ events generated using the 
same Monte Carlo event generator at the same y/s, one at hadron level, the other at detector level. 
Distributions normalised to the number of events at the detector and the hadron level are compared 
to derive bin-by-bin correction factors which are used to correct the observed x p distribution at each 
centre-of-mass energy. 

This bin-by-bin unfolding procedure is suitable for x p as the effects of finite resolution and accep- 
tance do not cause significant migration (and therefore correlation) between bins. Such a method is 
not readily applicable to multiplicity distributions, due to the large correlations between bins. Instead, 
a matrix correction is used to correct for detector resolution effects, followed by a bin-by-bin correction 
which accounts for the residual effects due to acceptance cuts and initial state radiation, as in previous 
OPAL multiplicity studies, e.g. [12,34]. 

Figure 3 shows the corrected values of mean charged particle multiplicity, (nt^), {n%f v ) and their 
difference A(n c h) as a function of centre-of-mass energy. It can be seen that, although the no-CR 
models vary in their predictions for both (n^) and (n^ they are in complete agreement that the 
value of A(n c h) is negligibly small, in contrast to the CR models shown in Figure 3(c). However, the 
predictions of both conventional QCD models and models of CR are found to be compatible with the 
data within uncertainties. As the multiplicity data are not observed to vary significantly with y/s, 
measurements from all centre-of-mass energies are combined assuming they are independent of y/s. 
The combined results are presented in Table 1, systematic uncertainties discussed in Section 4, and 
quantitative comparisons with models presented in Table 2. 
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3.2 Particle Flow 



The analysis of event properties presented here is a generalisation of the string effect analysis in 
three-jet e + e~ — > qqg events to the four-jet topology of W + W~ — ► qqqq. The situation is necessarily 
more complicated in the W + W~ — > qqqq channel because, in contrast to the three-jet case, events 
are not constrained by momentum conservation to be planar. The analysis is therefore carried out in 
four distinct planes, each of which is defined by a pair of jet axes. Charged particles and clusters of 
electromagnetic energy, selected as in [12], are combined into four jets using the k± [35] jet-finding 
algorithm, and the total momentum and energy of each of the jets are corrected empirically for 
double counting using the same energy flow algorithm employed in [20]. The jet momenta are further 
modified by a kinematic fit, imposing the four constraints of energy and momentum conservation, 
to obtain an improved estimate of the trajectories of the underlying four fermions from the W + W~ 
decays. To ensure that a relatively simple colour topology is being studied, events having a five-jet 
like topology 4 are rejected. This is also expected to lead to a better description of the dominant 
e + e~ — ► qq background events by the parton shower models. 

The analysis proceeds in three stages, namely: association of pairs of jets with W bosons and 
definition of four planes, projection of charged particles onto these planes, and comparison of the 
distributions of particles in these planes. Each of these aspects of the analysis is described below. 

The association between pairs of jets and W bosons is performed using a minor variant of the 
algorithm that was introduced in [27]. In the current scheme, the mass obtained from a five-constraint 5 
kinematic fit [20] is combined with the variables of [27] into a single likelihood discriminant, selecting 
the correct pairing of jets to W bosons with a purity of ~ 90%. The total number of events used 
in the analysis after all selections is 2199, with an overall efficiency for selecting W + W~ — > qqqq of 
~ 40%. Small variations in performance with centre-of-mass energy are detailed in Table 3. 

The pairing of jets originating from the same parent W defines two intra- W planes, as illustrated 
in Figure 4. With four jets, there exist two ways in which planes may be defined between jets which 
originate from different W boson (inter-W regions). The configuration which results in the smaller sum 
of inter- W angles is chosen. The motivation for this is the suggestion [2,4,5] that colour reconnection 
is more probable for topologies in which jets from different W bosons are close together in angle. In 
such a configuration, a rearrangement of the colour flow in the event would be energetically favoured 
due to a reduction in the overall "length" of the colour flux tubes. 

Reconstructed charged particles in the event are projected onto the intra- W and inter-W planes 
as follows, and illustrated in Figure 4. The first plane examined is that defined by the most energetic 
jet in the event ('jet 1') and the jet belonging to the same W ('jet 2'), as given by the jet pairing 
algorithm. The next plane considered is that between jet 2 and a jet ('jet 3') from the other W, such 
that the sum of inter-W angles is minimal. The third plane is the other intra-W region between jet 3 
and the remaining jet in the event, 'jet 4'. The final plane is that between jet 4 and jet 1. 

All charged particles in an event are projected onto each of the four planes in turn. In each plane, 
an azimuthal angle < xi < 2vr is defined, having positive sense between pairs of jets as described 
above and indicated in Figure 4. To account for the variation in the angle between pairs of jets, the 
distribution of particles is evaluated as a function of xi for each plane after rescaling, event-by-event, 
to the angle between the jets which define the plane, xo- This gives a normalised angle, xr = Xi/xo, 
where xr = corresponds to the jet axis of the lower number jet which defines the plane. Particles 

4 Following [20], five-jet events are classified as those in which the k± jet resolution parameter for the four-jet to 
five-jet transition, j/45, satisfies log(2/4s) > —5.6. 

The additional constraint imposed is equality of the masses of the two W boson candidates. 
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outside the inter-jet region, i.e., having \R > 1j are n °t considered further. In the case where a 
particle is projected into the inter-jet region of more than one plane, it is exclusively assigned to the 
plane relative to which it has the smallest transverse momentum. 

The four normalised inter-jet regions are combined in a single distribution in the range < % < 4, 
as shown in Figure 5(a), where the structure of the four jets is apparent. The variable % is defined as 
X = XR + ( n pianc — 1), where n p i anc is an integer between 1 and 4, corresponding to the four planes in 
the order given, and as shown in Figure 4. 

As seen in Figure 5(a), the data are consistent with the predictions of W + W~ production using the 
conventional hadronisation models plus the sum of all background processes. Figure 5(b) compares 
the sum of background-subtracted data with predictions from various CR models. The data are found 
to be adequately described by models, with the exception of the extreme scenario of the SK-I model, 
which predicts lower particle densities in the intra-W regions, and higher particle densities in the 
inter- W regions. 

In conventional QCD models without interaction of the colour fields between the W + and W _ , 
the particle density (or particle flow) is expected to be higher in the intra-W regions, < x < 1 an d 
2 < x < 3 than in the inter-W regions. After a rearrangement, in addition to a change in absolute 
number of charged particles in the event, there may be a migration of particle flow away from these 
regions in favour of the inter-W regions, 1 < x < 2 and 3 < x < 4. Consequently, one way of studying 
the effects of colour rearrangement is to compare the particle flow of the two intra-W regions to that 
of the two inter-W regions. As the properties of the two inter-W regions should both be affected by 
colour rearrangement in the same way, these are added together, as are the two intra-W regions. The 
ratio of intra-W to inter-W particle flow distributions is then formed, 



where n c h is the number of charged particle tracks projected into a given inter-jet region. 

Figure 6(a) compares the measured values of this ratio using all data (after background subtraction) 
with the predictions of conventional hadronisation models, while Figure 6(b) shows data compared 
with the predictions of various CR models. 

Differences between the data and models, and consistency between the predictions of different 
models, are more apparent in Figure 6 than in Figure 5. The data are found to lie slightly below the 
model predictions in the region away from the jet cores for conventional QCD Monte Carlo models 
and most CR models. While the SK-I model with kj = 100 shows significant separation from data and 
the other models, it is apparent that the predicted effects of CR in the other SK models are limited. 

3.2.1 Quantitative Measures of CR 

To quantify the consistency between data and predictions of models, the ratio of the integral of the 
particle flow in the intra-W regions to the integral of particle flow in the inter-W regions, 



is formed. This is a traditional observable used in string effect studies, e.g. [17], and is sensitive to 
differences in the number of particles in the inter-jet regions but relatively insensitive to their angular 
distribution and the choice of binning. 




& (intra-W) 

R 

^(inter-W) 



Rn 
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The limits of integration are chosen to optimise the predicted sensitivity in the SK-I model at 
y / s=189 GeV. This choice also allows the uncertainty on the ratio to be calculated from error prop- 
agation. In the case where the limits are extended too close to the cores of the jets, a significant 
correlation is introduced between neighbouring inter-jet regions and the error calculation is no longer 
valid. It is to be noted that to calculate the uncertainty correctly, the numerator and denominator 
of Equation 1 must be evaluated event-by-event, rather than by integration of distributions such as 
Figure 5. The validity of the statistical errors has been tested using data-sized samples from a variety 
of Monte Carlo models, each with more than 90 times the statistics of the entire 189-209 GeV data 
sample. 

To estimate the sensitivity to colour reconnection effects, the Rn predicted by each colour recon- 
nection model is compared with that obtained from the corresponding "no reconnection" scenario of 
the same model. As a guide to the performance of the analysis, the predicted statistical sensitivity of 
the analysis is summarised in Table 4 using model predictions at y/s = 199.5 GeV. The fraction of 
reconnected events in each model is also shown. The sensitivity is defined as the difference between 
a given reconnection model and its corresponding "no reconnection" sample (ARn), divided by the 
expected statistical uncertainty obtained using all data presented in this paper (cf^* ). It can be seen 
that there is a significant sensitivity to the extreme scenario of the SK-I model in which almost all 
events are reconnected but limited sensitivity to all other CR models considered. 

To combine the observed Rn from different centre-of-mass energies, an assumption has to be made 
about possible energy dependence of the measurements. Figure 7 shows the measured values of Rn 
together with the predictions of the JETSET and other Monte Carlo samples. Although some models 
exhibit an energy dependent Rn, the variations are small compared to the statistical uncertainties in 
data. Measurements are therefore combined assuming they are independent of y/s, and the impact of 
this is considered as a systematic effect. Quantitative comparisons of the combined Rn in data are 
made using predictions of all models studied at y/s = 199.5 GeV, at which energy a complete set of 
Monte Carlo samples is available. On the scale of variations predicted in Rn, this y/s is close to the 
luminosity-weighted centre-of-mass energy of 197.8 GeV. 



4 Systematic Uncertainties 

Systematic uncertainties are studied using measurements averaged over y/s, and are shown in Tables 1 
and 5 for the inclusive charged particle multiplicity and particle flow analyses, respectively. Sources 
of systematic error considered include hadronisation effects in the W + W~ models, detector effects 
related to tracking of charged particles and background subtraction. The analysis of mean particle 
multiplicities involves the unfolding of observed data to the hadron level and a possible additional 
uncertainty related to this procedure is studied. The particle flow analysis is performed using ratios 
of sums and no unfolding of the data is performed, and so many systematic effects are expected to 
cancel or be negligibly small. 

4.1 W+W - Hadronisation 

For the multiplicity analysis, samples of simulated W + W~ events incorporating JETSET hadroni- 
sation are treated as background-subtracted data and unfolded with each of the alternative W + W~ 
hadronisation models and colour reconnection models. The CR models are included in this procedure 
to allow for the possibility that such effects are present in the underlying physics. Where samples of 
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a given model exist at more than one ^fs the results of unfolding are averaged. The uncertainty is 
assigned as the largest difference between the values obtained when unfolding JETSET events with the 
default model (JETSET) and with any of the alternative models, and is dominated by the HERWIG 
models. 

For the particle flow analysis, this error is assigned by using samples of W + W~ events generated 
using KORALW, and hadronised with each of the models JETSET, HERWIG, ARIADNE, old JET- 
SET and AR-1. Where samples of a given model exist at more than one the results are averaged. 
Each model is treated as background-subtracted data, and the uncertainty is assigned as half of the 
maximum difference between the Rn predicted by any pair of models. This differs from the definition 
used for the multiplicity analysis: as no unfolding is performed, there is no default model against 
which to study systematic effects. It can be seen from Figure 7 that this uncertainty is determined by 
differences between the HERWIG and old JETSET hadronisation models. 

4.2 BEC 

While the presence of Bose-Einstein correlations among particles originating from the same hadroni- 
cally decaying W boson (intra-W BEC) has been unambiguously established [36], there is no significant 
evidence for BEC between particles originating from different W bosons (inter-W BEC) [37,38]. How- 
ever, these are not excluded and, following [38] and [20], a systematic error corresponding to 77% of 
the effect of (inter-W BEC) — (intra-W BEC) on the measurement is assigned. In the multiplicity 
analysis, the uncertainty was assigned as 77% of the difference in the hadron level multiplicity obtained 
when the inter-W BEC and intra-W BEC samples were each treated as background-subtracted data. 
In the particle flow analysis the uncertainty was assigned as 77% of the difference between the Rn 
values predicted by the inter-W BEC model and the intra-W BEC model. 

4.3 Track Definition 

Uncertainties arising from the selection of charged tracks are estimated by examining the stability of 
the difference between data and Monte Carlo predictions for multiplicity or Rn. Both analyses are 
repeated three times, with track selection criteria varied within reasonable limits [12]. The maximum 
allowed values of the distances of closest approach to the interaction region in r-(j> and z are varied 
from 2 cm to 5 cm and from 25 cm to 50 cm, respectively, and the minimum number of hits on tracks 
is varied from 20 to 40. The uncertainty on the charged track definition is the sum in quadrature 
of these three effects. This source represents a significant systematic effect and is dominated by the 
variation of the minimum number of hits required to form a track. 

4.4 Background 

Alternative models and cross-sections were used to estimate uncertainties associated with the back- 
ground subtraction. The uncertainty is formed using the difference between the measured multiplicities 
(or Rn values) obtained using the alternative background models and that obtained using the default 
background model and assumed cross-section at each centre-of-mass energy. 
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4.4.1 e+e — ► qq Modelling 

Uncertainties in generation of the hard process and hadronisation may affect the shape of the back- 
ground and are estimated by comparing models. A sample of e + e~ — > qq events, generated using 
KK2f and hadronised with each of JETSET, HERWIG and ARIADNE, is available at all centre-of- 
mass energies studied, as is a sample generated entirely using PYTHIA. This uncertainty is assigned 
as the largest difference between the result obtained using any model and the result obtained when 
the default JETSET model is used. 

4.4.2 e+e - — >■ qq Rate 

This uncertainty arises due to imperfect knowledge of the accepted background cross-section. It is 
evaluated using the deviations in the measurements caused when the e + e~ — > qq background rate is 
varied by ±5% and ±20% from its default value for W + W~ — > qqqq and W + W~ — > qq^£ events, 
respectively, where the allowed ranges are taken from [39]. 

4.4.3 Four-fermion Background Modelling 

This systematic uncertainty is estimated by using the KandY generator as an alternative to the 
default (KORALW) to simulate the WW-like four-fermion events, with the change in the measured 
multiplicity or Rjy value assigned as the uncertainty. Note that owing to the O{o) corrections in 
KandY, the alternative WW-like four-fermion cross-section is 2.5% lower than that of KORALW, and 
so the same 2.5% reduction is also applied to the cross-section of the KORALW W + W~ events when 
carrying out this test. 

4.4.4 Z°Z° Rate 

An uncertainty is assigned to the assumed cross-section for ZZ events. It is estimated by varying the 
ZZ component of the four-fermion background by ±11% for W + W~ — > qqqq events [40], and ±20% 
for W + W _ — ► qqft^ events [39]. 

4.4.5 Residual Backgrounds 

Two further small sources of background are considered. The first small source of background is only 
relevant for W + W~ — ► qqiT>e events and so is considered for the multiplicity analysis alone. It is 
assigned as the effect observed on the measurements when all predicted e + e~ — > l + t~ backgrounds 
are neglected. 

The second source is due to two-photon background, and is estimated as the difference found in the 
final result when the small, predicted background from this source is included (default) or neglected. 
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4.5 Unfolding Method 



For the multiplicity analysis, the results of the x p and direct multiplicity analyses are compared and 
the difference in their central values is assigned as a source of possible uncertainty. 

4.6 Centre-of-Mass Energy Dependence 

The multiplicity and particle flow measurements are assumed to be independent of y/s, as discussed 
in Section 3. An alternative choice considered for the particle flow analysis was to correct the mea- 
surements of Rn according to the weak energy dependence predicted by the KORALW model with 
JETSET hadronisation. The difference between these two assumptions is found to be small, at a level 
of 2% of the statistical uncertainty on the combined result, and is therefore neglected. 

4.7 Cross-Check Using W+W~ — > qq£T^ Data 

By way of a cross-check that the data are adequately described by the conventional hadronisation 
models, the particle flow analysis is repeated using W + W~ — ► qqiue events, in which there can be no 
(inter-W) colour reconnection. The event selection is restricted to events in which the charged lepton 
is classified as either an electron or a muon. The four planes used in these events are defined by the jet 
or fermion directions derived from a kinematic fit in which the constraints of energy and momentum 
conservation are imposed (4-C fit, as in [20]). Figure 8(a) shows the distribution of particle flow for 
all W + W~ — > qqtPe data, which are found to be described well by the predictions of the JETSET, 
HERWIG, ARIADNE and old JETSET models. 

By construction, the two jet axes corresponding to the hadronically decaying W boson are centred 
at x = and % = 1, while the direction of the charged lepton and that inferred for the unobserved 
neutrino are at \ = 2 and % = 3, respectively. Note that charged particles associated with the 
leptonically decaying W bosons are not included in these figures. The non-zero multiplicity in the 
region between the two leptonic "jets" is due to the particles projected into this plane from the 
hadronically decaying W. Similarly, the abrupt change in the distribution in the region of the leptonic 
W is because none of the charged particles projected onto this plane plays a direct role in defining it. 

Figure 8(b) shows the energy evolution of Rn for W + W~ — ► qqiT>e events, together with the 
predictions of the same set of four hadronisation models. Similarly, the models provide a reason- 
able description of the data, giving confidence to the analysis procedure. No additional systematic 
uncertainty is assigned as a result of this study. 



5 Results and Discussion 

The measurements of the mean charged particle multiplicities corrected to the hadron level and aver- 
aged over the range y/s ~ 189-209 GeV, are: 




38.74 ±0.12 ±0.26 



19.39 ±0.11 ±0.09 



A(n ch ) = 



-0.04 ±0.25 ±0.17 , 
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where in each case the first uncertainty is statistical and the second systematic. The difference in mean 
charged particle multiplicities in hadronic W decays in qqqq and qqtp£ events, A(rj c h), is found to be 
consistent with zero within uncertainties. All models considered are found to lie within 1.2 standard 
deviations of the measurement, as shown in Table 2. As no evidence is found in this measurement 
for colour reconnection between the two hadronically decaying W bosons, the average of data from 
W + W _ — > qqqq and W + W _ — > qqft^ events, weighted by statistical uncertainties and taking into 
account correlations in the systematic uncertainties, is used to yield a measurement of the charged 
particle multiplicity from a single hadronically decaying W, 

(n^) = 19.38 ± 0.05(stat.) ± 0.08(syst.) . 

As this average is made under the assumption that there is no colour reconnection between W bosons, 
the CR contribution to the W + W~ hadronisation uncertainty has been removed. 

The analysis of particle flow allows a simple comparison with models of colour reconnection, using 
the data of Table 5. The measurement obtained using approximately 625 pb _1 of data in the range 
189-209 GeV yields: 

R N = 1-243 ± 0.025(stat.) ± 0.023(syst.) . (2) 

This result may be compared with the predictions of the models at y / i=199.5 GeV given in Table 2. It 
can be seen that Rn measured in data is lower than all models except the SK-I (kj = 100) sample. The 
(signed) significance of these differences is also presented, varying from approximately 4.4 standard 
deviations of the total error (cr^ al ) for an extreme scenario of the SK-I model, to -2.0cr^* al for 
HERWIG, with most other models populating a region around — lcr^ al . 

Comparing the measured Rjy with the predictions of Table 2, the data are seen to be closest to 
the predictions of the SK-I model with strength parameter kj = 0.9. As this parameter is arbitrary, 
it may be varied to optimise the consistency with the measured Rn of Equation 2. The Ax 2 curve 
corresponding to this variation is presented in Figure 9. Parametrising this curve using a fourth 
order polynomial, the best agreement with data is obtained when approximately 37% of events are 
reconnected in the SK-I model, corresponding to the value kj = 1.0. The 68% confidence level allowed 
region deduced from the A% 2 curve corresponds to 0.10 < p rcco < 0.56. This result is not combined 
with the analysis of inclusive particle multiplicity as they have correlated systematic uncertainties and 
no significant improvement in sensitivity is expected. 

It should be noted that the properties of the SK-I model vary significantly with the parton shower 
cut-off parameter and therefore this kj cannot be directly compared to similar results from other LEP 
collaborations. Any combination of results from the different experiments is best performed on the 
basis of analysis of a common set of simulated events, analysed independently by each experimental 
collaboration [18]. 



6 Conclusions 



The predictions of models of colour reconnection implemented within the ARIADNE Monte Carlo, 
the HERWIG model and the SK model, have been compared with OPAL data recorded at ^/s ~189- 
209 GeV using both inclusive measurements of particle multiplicity and a generalisation of the "string 
effect" analysis to the four-jet topology of W + W~ — > qqqq events. 

Studies of reconnection phenomena using the extreme scenarios of the SK-I model show that 
changes up to approximately 1% may be expected in (nJJ), where the total experimental uncertainty 
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on measurements of (n ?) is 0.7%. Other models predict somewhat smaller effects. Defining A(n c h) 
using data alone provides a model-independent (but less sensitive) test of possible reconnection effects. 
The inclusive measurements of particle multiplicity find no evidence for such effects. 

Measurements of particle flow in the OPAL data exclude an extreme scenario of the SK-I model 
and are compatible with other CR models such as SK-II, SK-II', AR-2 and that of HERWIG. They 
are also compatible with models which do not include colour reconnection, slightly disfavouring the 
conventional HERWIG model. The results of this analysis are not combined with measurements 
of inclusive particle multiplicity as they have correlated systematic uncertainties and no significant 
improvement in sensitivity is expected. The best agreement with data is obtained using the SK-I model 
with a reconnection probability, p rcco ) of approximately 37%, corresponding to a model parameter ki = 
1.0 within the context of the OPAL tuning of the JETSET hadronisation model. The 68% confidence 
level allowed region deduced from the x 2 curve corresponds to 0.10 < P rcco < 0.56. This result is used 
to help constrain the systematic uncertainty related to colour reconnection in measurements of the W 
boson mass [20]. 
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— U.U4 
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ly.oo 


Stat, error 


n 1 9 

U. 1Z 


nil 


n 9^ 


n n^ 


Systematics 










W + W~ hadronisation 


0.22 


n no 

0.08 


O.Oo 


O.Ofa 




n no 

0.03 


0.01 


0.05 


n m 

0.01 


Irack definition 


n no 

u.uy 


O.Oo 


n nn 

0.09 


0.04 


e^e — ► qq modelling 


nil 


n m 


n 1 n 


n n/i 




n ni 

VJ.VJ ± 


n nn 




n ni 

VJ.VJ ± 


Four-fermion background modelling 


0.01 


0.02 


0.02 


0.01 


Z°Z° rate 


0.01 


0.01 


0.02 


0.00 


Residual backgrounds 


0.00 


0.01 


0.02 


0.00 


Unfolding procedure 


0.00 


0.01 


0.02 


0.00 


Total syst. 


0.26 


0.09 


0.17 


0.08 



Table 1: Results and estimated systematic effects in inclusive charged particle multiplicity measure- 
ments, see text for details. 



Sample 


A(n ch ) 


significance (A (n c h ) ) 


Rn 


significance(i?jv) 


Data 


-0.04 ±0.30 




1.243 ± 0.034 




SK-I(fc/ = 100) 


-0.42 


+1.2 


1.092 


+4.4 


SK-I(fcj = 0.9) 


-0.29 


+0.8 


1.246 


-0.1 


SK-II 


-0.14 


+0.3 


1.273 


-0.9 


SK-II' 


-0.16 


+0.4 


1.277 


-1.0 


AR-2 


-0.19 


+0.5 


1.271 


-0.8 


HERWIG-CR 


+0.32 


-1.2 


1.282 


-1.2 


JETSET 


-0.04 


0.0 


1.291 


-1.4 


HERWIG 


+0.02 


-0.2 


1.311 


-2.0 


ARIADNE 


+0.02 


-0.2 


1.286 


-1.3 


Old JETSET 


-0.03 


0.0 


1.280 


-1.1 


AR-1 


0.00 


-0.2 


1.304 


-1.8 



Table 2: Comparison of the average measured A(n c h) and Rn in data with various models at 
\fs = 199.5 GeV. The level of agreement in Rn is given by the significance (R n), i-e., (Rn {data) — 
R n (model) )/c^ r aI , the difference between the average value of Rn in data and each model divided 
by the total uncertainty, and similarly for A(n c h). 



(y/s) (GeV) 


jCdt (pb- 1 ) 


Selected events 


Efficiency (%) 


Purity (%) 


Correct jet pairing (%) 


188.6 


183.0 


675 


39.7 


86.1 


90.3 


191.6 


29.3 


92 


39.1 


86.5 


90.0 


195.5 


76.4 


277 


40.0 


87.1 


89.9 


199.5 


76.6 


253 


38.9 


87.0 


89.4 


201.6 


37.7 


145 


38.4 


84.2 


89.1 


206.0 


220.5 


757 


37.7 


86.2 


88.5 



Table 3: Summary of the integrated luminosity and number of candidate events used in the particle 
flow analysis, after all selection criteria, at each centre-of-mass energy in the range 189-209 GeV. The 
efficiency and purity are defined relative to the W + W~ — ► qqqq production process. The "correct" 
pairing in the rightmost column is defined by whichever association of observed jets in the detector 
minimises the sum of angular differences relative to the original four fermions from the W + W~ decay. 
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Model 


prcco ( % ) 


AR N /at:~ 


SK-I(fc/ = 100) 


98.2 


7.9 


SK-I(fcj = 0.9) 


34.3 


1.7 


SK-II 


17.2 


0.6 


sk-ii' 


16.1 


0.5 


AR-2 


49.4 


1.3 


HERWIG-CR 


23.0 


0.9 



Table 4: Summary of the predicted statistical sensitivity of the particle flow analysis for different 
models of colour reconnection at y/s = 199.5 GeV. The sensitivity is defined as the difference between 
a given reconnection model and its corresponding "no reconnection" sample (ARn), divided by the 
expected error of all data combined (cfj^ - )- P rcco is the fraction of colour reconnected events in each 
model. For AR-2, the no-CR model is AR-1. 



Particle Flow 


Rn 


Data 


1.243 


Stat, error 


0.025 


Systematics 




W + W~ hadronisation 


0.015 


BEC 


0.002 


Track definition 


0.014 


e + e~ — ► qq modelling 


0.010 


e + e~ — ► qq rate 


0.002 


Four-fermion background modelling 


0.002 


Z°Z° rate 


0.001 


Residual backgrounds 


0.000 


Total syst. 


0.023 



Table 5: Result and estimated systematic effects in particle flow measurements, see text for details. 
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Figure 1: Uncorrected charged particle multiplicity distributions for data in the range ^/s=189- 
209 GeV: (a) W + W~ — > qqqq events and (b) the hadronic part of W + W~ — > qq^£ events. Points 
indicate the data with statistical errors, lines show the expected sum of signal and background con- 
tributions for a variety of signal models, and the hatched histogram shows the expected background. 
Predictions of the conventional QCD hadronisation models JETSET and HERWIG, the AR-1 model 
and the 100% CR SK-I model, are shown. 
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Figure 2: Uncorrected x p distributions for data in the range -^75=189-209 GeV: (a) W + W~ —>■ qqqq 
events and (b) the hadronic part of W + W~ — > qq^v events. Points indicate the data with statistical 
errors, smooth curves show the expected sum of signal and background contributions for a variety of 
signal models, and the hatched region shows the expected background. Predictions of the conventional 
QCD hadronisation models JETSET and HERWIG, the AR-1 model and the 100% CR SK-I model, are 
shown. Monte Carlo samples are normalised to the predicted number of signal plus background events, 
therefore the hatched regions correspond to the mean number of particles in candidate W + W~ events 
which originate from background sources, rather than the mean number of particles per background 
event. 
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Figure 3: Centre-of-mass energy dependence of the measured (unfolded) mean charged particle 
multiplicity for (a) W + W~ — > qqqq events, (b) W + W~ — ► qqft^ events, and (c) the difference, 
A(n c h) = («ch) — ^( n ch £U )- Points indicate the data with statistical errors and lines the predictions 
of W + W~ models incorporating either conventional QCD hadronisation or CR. The predictions of 
JETSET, HERWIG, ARIADNE, the old tune of JETSET and AR-1 are indistinguishable from zero 
in (c), in all cases having values smaller in magnitude than 0.05, and so are not shown. 
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Figure 4: Illustration of inter-W and intra-W regions, and numbering convention adopted, as described 
in Section 3. 
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Figure 5: The particle flow distribution in the four inter-jet planes, normalised event-by-event to the 
inter-jet angles, as described in the text and in Figure 4. Points represent the data with statistical er- 
rors, (a) compares the data with the predictions of conventional QCD hadronisation Monte Carlo mod- 
els and AR-1, (b) compares the data, after background subtraction, with several CR models. Monte 
Carlo samples are normalised to the predicted number of signal plus background events, therefore the 
hatched region in (a) corresponds to the mean number of particles in candidate W + W~ — > qqqq events 
which originate from background sources, rather than the mean number of particles per background 
event. 
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Figure 6: Distribution of i?fl OW ) the ratio of the particle flow in the intra-W regions to that in the 
inter-W regions. Points represent the data after background subtraction, with statistical errors, (a) 
compares the data with the predictions of conventional QCD hadronisation Monte Carlo models and 
AR-1, (b) compares the data with CR models. Note that there are correlations between bins in these 
distributions. 
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Figure 7: Centre-of-mass energy dependence of Rn for W + W~ — > qqqq events. Points indicate the 
data with statistical errors and lines the predictions of W + W~ models incorporating either conven- 
tional QCD hadronisation or CR. The predictions of the HERWIG CR model, AR-1 and AR-2, which 
are only available at -^=199.5 GeV, are indicated by symbols. 
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Figure 8: (a) The particle flow distribution (as in Figure 5) for W + W~ — ► qqiUg events and con- 
ventional QCD hadronisation Monte Carlo models, as described in Section 4.7. Points represent the 
data with statistical errors and the hatched region the sum of all background contributions, (b) The 
energy evolution of Rn as measured in W + W~ — ► qqiVi events, in comparison with the predictions 
of conventional QCD models. 
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Figure 9: Ax 2 curve obtained from comparison of the average Rn measured using OPAL data between 
189 GeV and 209 GeV, with the predictions of the SK-I model as a function of the fraction of 
reconnected events, p Teco i carried out at a reference centre-of-mass energy of 199.5 GeV. (a) shows 
the entire range of _p rcco ) while (b) shows the lower p reco range of (a) in more detail. The best 
agreement between the model and data is obtained when 37% of events are reconnected in this model. 
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